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We reconsider the problem of the spontaneous emission of light by an excited atomic state. We 
scrutinize the survival probability of this excited state for very short times, in the so-called Zeno 
regime, for which we show that the dynamics is dictated by a coherent -in phase - response of the 
on-shell and off-shell vacuum modes. We also develop a perturbative approach in order to interpolate 
between different temporal regimes: the Zeno, golden rule (linear) and Wigner-Weisskopf (exponen¬ 
tial) regimes. We compare results obtained with the E • x and A ■ p interaction Hamiltonian,s, using 
successively the dipole approximation and the exact coupling. 
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I. INTRODUCTION 

We study the spontaneous emission of light by an 
atomic dipole (two-level quantum system) at zero tem¬ 
perature in which an electronic state is initially excited 
while all the electromagnetic modes are taken to be in 
their vacuum (zero-photon) state. In the standard treat¬ 
ment of this kind of problem, one derives the Fermi 
golden rule jT|, which states that the probability that 
the dipole stays in its excited state (the so-called sur¬ 
vival probability) decreases linearly with time. At longer 
times, the standard treatment is the Wigner-Weisskopf 
approach [2], in the framework of which the survival 
probability is predicted to decrease exponentially. It is 
well known [3], on the other hand, that, under very gen¬ 
eral circumstances, the survival probability of a quantum 
system decays quadratically -and not linearly- with time 
at very short times. This initial quadratic behaviour is 
known as the quantum Zeno regime. 

We studied in detail the behaviour of an atomic dipole 
in the Zeno regime. Our calculations predict a small 
but sudden decrease of the survival probability, in con¬ 
junction with the emission of off-resonant light. For the 
2p — Is transition in atomic hydrogen -within the rotat¬ 
ing wave approximation- we found that electromagnetic 
modes with frequency up to a thousand times larger than 
that of the atomic transition are excited. 

As we shall show, the transient Zeno regime is dom¬ 
inated by the collective, in phase, response of the elec¬ 
tromagnetic field modes. Very quickly, this collective 
response gets out of phase (incoherent) so that non¬ 
resonant electromagnetic modes become, in a way, in¬ 
active. It is only at this point that the linear decay of 
the survival probability arises, as described by the Fermi 
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golden rule. This linear decay is essentially driven by the 
resonant modes of the electromagnetic field. 

Our paper is organized as follows. In Sect. [TT] we sum¬ 
marize the standard quantum electrodynamical (QED) 
description of dipolar decay, including the Fermi golden 
rule. In Sect. |III| we describe a perturbative approach 
to the problem in which the outgoing electromagnetic 
modes are discretized. In this approach only nearly res¬ 
onant (on shell) modes are considered. This approach 
makes it possible to simulate an exponential decay of the 
survival probability in full agreement with the Wigner- 
Weisskopf approximation. It also predicts a Zeno effect 
which can be considered to be “standard”: for short times 
the parabolic dependence in time of the survival prob¬ 
ability results in a short “plateau” directly followed by 
a linear time dependence (Fermi golden rule), naturally 
merging at longer times into the exponential behaviour. 

In Sects. [TVl |V| and |Vl| we generalize the perturbative 
approach of Sect. m by taking account of the response 
of off-resonant electromagnetic modes. For very short 
times, their collective response results in a non-standard 
Zeno effect which is the essential new result of our paper. 
We first study this ef fect in the case for the E • x atom- 
field coupling (Sects. IV and |~vj), and then for the A ■ p 
coupling (Sect. VI). In both cases, the collective answer 
of the off-resonant modes induces a small decrease of the 
survival probability for short times, leading to the usual 
Fermi and Wigner-Weisskopf regimes. The last section is 
devoted to open questions and conclusions. 


II. BASIC ELEMENTS ON THE DECAY OF A 
QUANTUM DIPOLE 

A. Hamiltonian of the atom-light interaction 

We consider a two-level atom (ground state | g), excited 
state |e)) interacting with the electromagnetic field. The 
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atom sits either in free space or in an ideal metallic cav¬ 
ity. In both cases, the field can be expanded in terms of 
normal modes labelled by an index A. The Hamiltonian 
reads 


H = H a + H r + H,. (1) 

where 

Ha = E s | g) (g | ~\~F e | e) (e | (2) 

is the atomic Hamiltonian ( E g and E e being the energies 
of the ground and excited state), 

Hr = A ^d A + 0 (3) 

is the Hamiltonian of the free radiation (a A is the anni¬ 
hilation operator for photons in the mode A and w A the 
angular frequency of these photons). Let us consider for 
the moment that the atom-field interaction Hamiltonian 
reads 


H t = - d • 


DM 

eo 


( 4 ) 


which is valid in the electric dipole approximation. Here, 
d is the atomic dipole operator and D (r 0 ) is the dis¬ 
placement field at the atom location ro- We shall go 
beyond this approximation later but for the time being 
we prefer to formulate the problem at this level of ap¬ 
proximation, having in mind that the techniques that we 
shall develop in the next sections can easily be general¬ 
ized to other choices of interaction Hamiltonians. Within 
the two-level atom approximation, the dipole d can be 
written as 


d = F d il = I Ml d \j){j\= X]dy U)0' I. ( 5 ) 

i,j i,j 


where (i,j) £ {e,g} 2 and dy = (* | d | j). Because 
of spherical symmetry, d gg = d ee = 0 and d ge = d* g . 
Therefore the electric dipole operator can be written 


d = d ge |g)(e| +d* e |e)(g| . (6) 

Since there are no free charges in the system, the dis¬ 
placement held is purely transverse everywhere in space 
(in contrast to the electric held which is not transverse at 
ro) and can be expanded on the transverse modal func¬ 
tions v R as 


D(r) 

eo 


1 51 (® A Va “ 4 V I* ( r )) ( 7 ) 


if the modal functions are normalized such that 


J v[(r)’V?(r) dr = S\y. 


( 8 ) 


By plugging the expansions (|6| and ([TJ) into Q, we can 
write the interaction Hamiltonian Hi as the sum of two 
contributions H R and Hf R , where 

H? = ihJ2(G\( r o) 4 |g)( e l ~ G \ ( r o) a x |e>(g|) 

(9) 

contains the rotating interaction terms and 

Hf R = ihJ2(F x ( r 0 ) 4 | e) (g | —-F A (r 0 ) a A |g)(e|) 

( 10 ) 

contains the counter-rotating interaction terms. The 
coupling constants G\ (r 0 ) and F\ (r 0 ) are given by 


(ro) = dge ■ v£* (r 0 ), (Ha) 

Fx (r 0 ) = \j^~ d ge ' v a* ( r o) • (lib) 


The interaction terms appearing in H R M are easy to 
interpret. The hrst term corresponds to the creation of 
a photon while the atom jumps from its excited to its 
ground state. The second term represents the reverse 
process, that is, photon absorption. 

The interaction Hamiltonian Hf R (101 is sometimes 
claimed not to conserve energy [3]. The hrst term rep¬ 
resents the process of adding one photon to the held while 
simultaneously promoting the atom from the ground to 
the excited state; the second one represents the reverse 
process. In the derivation of the Fermi golden rule as 
well as in the Wigner-Weisskopf approach, it is common 
to neglect the counter-rotating interaction terms Hf R , 
which is called the rotating-wave approximation (RWA) 
M- As our goal is ultimately to estimate corrections 
to the W-W predictions, we shall from now on conhne 
ourselves to the same regime and thus stick to the Ro¬ 
tating Wave Approximation. This is particularly justi¬ 
fied in studies of the short-time behaviour of the system, 
because to hrst order in perturbation theory, the counter¬ 
rotating interaction Hamiltonian (10) does not contrib¬ 
ute. On the other hand, we shall not neglect the emission 
of off-resonant photons as is traditionally done. 


B. Decay equations within the Rotating Wave 
Approximation 

We now consider the problem of an atom that is ini¬ 
tially in its excited state | e) and decays to its ground 
state | g) due to its coupling to the electromagnetic held. 
If the held is initially unpopulated, the quantum state of 
the system at any time can be written as 

\ip (t)) = c e (f) e _ s(^+ E vac)t |e, 0) 

+ 5>sM) e-Me+^+Wac)* | g , 1a ), ( 12 ) 

A 
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where | e, 0) means that the atom is in the excited state 
and the field contains no photons and | g, 1;\) means that 
the atom is in the ground state and the field contains a 
single photon in the mode A. In the absence of any inter¬ 
action (Hi = 0), the coefficients c e and c gj \ do not evolve 
and each term of the superposition oscillates at its own ei- 
genfrequency. In the so-called weak coupling regime, the 
functions c e and c gt \ are expected to be slowly varying 
in time. The vacuum energy E vac = hu>\/2 produces 

a phase shift oscillating at an infinite angular frequency 
E vac /h. However, this phase shift can be factored out and 
appears as a global phase that can be simply ignored. 

By injecting (12) in the Schrodinger equation 


iH dt = {Ha + H r 


H?) | 


(13) 


we find the differential equations governing the evolution 
of functions c e and c gi A- They read 


c e = M c gjA (t) e 1( “* Weg)t , (14a) 

A 

c g ,A (t) = Ga (r 0 ) Ce (t) (14b) 

Up to the RWA, these equations are exact and apply to 
both emission in free-space and cavity emission. They 
possess no known exact solution pQ. 

Wigner and Weisskopf [2] attempted to solve these 
equations. They analysed the case of spontaneous emis¬ 
sion in free space using the ansatz 

c e (f)«e - S‘ (15) 


to explain the universally observed exponential decay of 
the excited state population. This solution is not an ex¬ 
act solution, but is widely considered to be a satisfactory 
approximation at long times. Using this ansatz, Wigner 
and Weisskopf found that T is well approximated by 


r = 


W eg l^gel 

3nheoC 3 


(16) 


which coincides with the value found by applying the 
Fermi golden rule, as we now discuss. 


C. Fermi golden rule from time-independent 
perturbation theory 

Consider the interaction Hamiltonian Hj as a per¬ 
turbation. The Fermi golden rule is usually obtained 
from time-dependent perturbation theory but here we 
show how time-independent perturbation theory to the 
first order also allows to retrieve the Fermi golden rule. 
This original approach opens the way to a more general 
time-independent perturbative treatment of the problem 
that will be developed later. The eigenstates of the free 
Hamiltonian Ha + H r are | e, 0) and | g, 1a), with respect¬ 
ive eigenenergies E e = hui e and E g + hui\ = h (co g + ui\). 


Perturbation theory pj tells us that, treating Hj to first 
order, these eigenenergies are unmodified. The normal¬ 
ized perturbed eigenstates are 


5 , 0 )<» = -b(|e, 0 ) + i5:Ad£sL| &1 A 
VN y “ w eg - WA J 


|g,U> (1) = 


V^Va 


g,lx)+i Gl (ro) |e, 0) 


^eg 


where 


N = l + J2 


IGaMI 2 

(W eg ~ Wa) 2 


N x = 1 + 


|Ga (r 0 )P 

(W eg - Wa) 


2 ’ 


(17a) 

(17b) 

(18a) 

(18b) 


Note that the perturbed eigenstates © are orthogonal 
to lowest order in the perturbation. Notice that, to lowest 
order, we have 

M> = ^ (|e,0><‘> -i|g, U >(»' 


N uj eg uj\ 


(19) 

We can thus write the evolution of the system when the 
atom is initially excited, with no photon in the field. Us¬ 
ing the spectral theorem for the total Hamiltonian H 
with its approximate eigenstates © derived from time- 
independent perturbation theory, we get 


lUJpg-t 


(e,0|e-^‘ |e,0) =- 


N 


,2 e 1 


i {uie S -Ul X )t 


1 + ^2 i g a ( r o)r. 2 

v A K g -W A ) , 


( 20 ) 


After a bit of algebra, this yields 
2 


(e, 01 e~^ Ht I e. 




N 2 ^ (w eg - W A ) 2 


sin 2 (w eg - w A ) 


t 
2 
( 21 ) 

to lowest order. Let us define Pdecay as 1 minus the sur- 

i - 2 

vival probability: Pdecay (t) = 1 — (e, 01 e~^ Ht | e, 0) . 
In free space, the sum over free modes is an integral. 
The substitution follows the procedure 




X=1 * 


d 3 fc 

(2^ 


( 22 ) 


where x are the two polarizations, V is the initial quant¬ 
ization volume and the modal functions are given by 


(r) —> (k, r) = 


Q ik-r 




Hx) ( k ) • 


(23) 
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The unit polarization vectors obey [7] the closure relation 

k l ki 


2 

E 

>c= 1 


(E)) (k)e^(k)=^- 


k 2 


(24) 


With all this in mind we compute the decay probability 


P, 


decay 


(t) = 


4 Vc 

2hen 


f d 3 k 

dge ' £ (x) (k) 

/ (27t) 3 

V 


X Ikl 


sin 2 (Kg ~ c[k|) |) 
(w eg -c|k|) 2 


2c 

f d 3 k 

Id I 2 

d • k 

2" 

free J 

(2tt) 3 

l a ge| 

ge |k| 



X Ikl 


sin 2 ((w eg - cjk|) 


t 1 


g e I 


67r 2 eo he 3 


(W eg — C |k|) 

/*oo 

/ do; w 3 


sme 


(w eg w) 


(25) 


where we performed the angular integration by requir¬ 
ing that the polar angle 9 be zero when k points in the 
direction of d ge . 

In the limit of large times, that is, times larger than 
l/w e g, one finds the usual Fermi golden rule with the 
decay rate T (16), that is, Pdecay (t) = Ft. In this regime, 
the cardinal sine in (25) becomes very peaked around 
the atomic transition frequency. One can then consider 
that only frequencies close to w eg contribute to (251 and 
rewrite the integral in the latter equation as 


/ 

J UJ 


tdeg+A 


do; w 3 sine 2 


„ —A 


(W eg - U) ^ 


Then, making use of 


/*CD e g +A 

(W eg - W) l 

/ dw w 3 sine 2 

Ju><, g -A 

v to 2 


p + OO 


t —^ -(“00 


= UJZ 


W eg 

27r 

T’ 


d ui sine 2 


(w eg w) 


(26) 


the Fermi Golden Rule is established. This is only valid 
at so-called intermediate times, larger than l/w eg , but 
smaller than 1/T, so that the perturbative approach still 
makes sense. 

Strictly speaking, taking the limit t —> +00 under the 
integral on the right-hand side of (25), as is usually done 
to obtain the Fermi golden rule, is not rigorous, since 
that integral diverges (for all times). For instance, for 
sufficiently small times the cardinal sine in (25) can be 
taken to be equal to 1, and hence the integral behaves 
like J Q +0O dww 3 = + 00 . To accommodate this, we in¬ 
troduce a cutoff on frequencies. A crude but natural [5; 


cutoff frequency is the one corresponding to the size of 
the atomic dipole we considered: Kq = 2e7r/ |d ge | with e 
the elementary charge. Note that in a recent treatment 
of a similar problem |2], the Compton frequency was used 
as a cutoff. We shall see later that when one goes beyond 
the dipole approximation, a natural cutoff appears from 
the computations and need not be implemented artifi¬ 
cially (it is of the same order of magnitude as A'c). We 
then have, for very small times, 


P 


decay 


(t) = f2 |dgc 


rcK c 


67r 2 eo he 3 
c 


dwe 


IQ 


24^.* ' d -' 2 f2 - 


(27) 


A parabolic Zeno behaviour is thusly obtained. To con¬ 
nect this initial behaviour with that, at longer times, 
predicted by the Fermi golden rule and the Wigner- 
Weisskopf approximation, we make appeal to numerical 
calculations. As a prerequisite, we now investigate how 
results from Wigner and Weisskopf are retrieved, in a dis¬ 
cretized approach in which the two-level atom is coupled 
with a finite number of quasiresonant modes of the elec¬ 
tromagnetic field. 


III. EXPONENTIAL DECAY FOR A 
TWO-LEVEL ATOM COUPLED WITH 
QUASIRESONANT FIELD MODES 


A. A discrete approach to spontaneous emission 


Instead of a continuum of field modes, we will now 
consider a finite amount thereof. We will thus pass from 
a continuous regime where Fermi’s Golden rule can be 
applied to a discrete regime where the sum in (12) is 
both discrete and finite. The Hilbert space is accordingly 
assumed to be finite-dimensional, which opens the way 
to a considerably simplified numerical treatment based 
on matrix equations. 


The problem that we face is twofold: for each energy 
there is a continuous infinity of field modes which cor¬ 
responds to all possible spatial directions of emission of 
a photon; moreover we must consider a range of ener¬ 
gies running from 0 to HcKq. Usually, only resonant 
(on-shell) photons are considered, and the multiplicity 
of outgoing photonic modes is hidden in the Density of 
States (DoS) formulation of the Fermi golden rule Jj. We 
shall adopt here another strategy, in which we integrate 
over the continuum of spatial directions of emission (and 
over various polarizations) without resorting to any kind 
of approximation. 


In order to explain the principle of our method, let us 
first consider a situation in which two electromagnetic 
modes of equal energy only are present. The total (inter- 
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action + free) and free Hamiltonians then read 



huio 

-i hG\ 

-ihGi 


huio 

0 

0 ' 

H = 

ifiGi 

hu 

0 

H 0 = 

0 

hu 

0 


}hG 2 

0 

hu 


0 

0 

hu 


The corresponding Schrodinger equation is 


Xo (t) 


x 0 (t) 

X\ (t) 

= H 0 

Xi (t) 

x 2 (t) 


X 2 {t) 


+ i h 


-G\x i (i) - G* 2 x 2 (t) 

Gix 0 (t) 

G 2 x 0 (t) 


(29) 


Let us now define the states |A) and | B) through 


W 
I B) 


i 

V\Gi\ 2 + \G 2 \2 
i 

V\Gi\ 2 + \G 2 \2 


0 

G\ 

G 2_ 

0 

g* 2 

-G\ 


(30) 

(31) 


In the new orthonormal basis spanned by the states 
(1 0 0) T , | A) and | B), the Schrodinger equation now reads 


X 0 (t) 


Geff X A ( t ) + W 0 X 0 (t) 

XA (t) 

= 

Geff X 0 (t) + UX A (t) 

XB (t) 


UX B (t) 


(32) 


where w e introduced the effective coupling constant 
G e ff = \J\Gi\ 2 + |G 2 | 2 . Obviously, the interaction only 
involves the state |A), while the orthogonal state | B) ig¬ 
nores the coupling with the “excited atom” state (1 0 0) T . 
This is still true even if an arbitrary number of modes 
(say N) are coupled to the excited state: the states or¬ 
thogonal to 


\A) 


1 

V\Gi\ 2 + \G 2 \ 2 + --- + \G n \ 2 


- 0 ■ 
G\ 
G 2 


G n . 


(33) 


numerically integrate over slices of various energies, we 
have to discretize the energy, in order to implement an 
approximation scheme that we can tackle numerically. 

Making use of equations ©, (|nj, @ and @, it is 
easy to check that a slice of energy comprised between 
hck\ and hck x+ \ is characterized by an effective coup¬ 
ling constant equal to 


^cff — 


2 /*27T /■ 7T 

' dip I d0 sin0x 
o Jo 


E 

.>c=l 


rfc " +1 k 2 dk 2 

—3 |G„ (k, r 0 )| 

fcA (2 tt) 


(35) 


As shown in App. 1 2 a| where the effective Hamiltonian is 
rewritten in terms of (averaged) ladder operators, if we 
successively perform the integration over the angles and 
over energies, we find 




(36) 


We also fixed (see App. 2a) the energy of the ef¬ 


fective mode assigned to an energy slice by impos¬ 
ing that it be equal to the average energy over the 
slice, weighted by the DoS. This yields, for an en¬ 
ergy slice bounded by k\ and fc A+ i, the effective energy 
= (4/5) he (k 5 x+1 - k \) / (k{ +1 - kf ). From now on, 
we drop the e ff subscripts from these quantities. The 
Hamiltonian of the system is as before of the form 


H = hu eg |e, 0><e,0|+£ hu x |g, 1 A ) (g, 1 A | 

A 

+ i^E [ Ga ( r °) ls» °l “ g a ( r o) |e, 0)(g, 1 A |], 

(3 7 ) 

excepted that here again summations run over a finite 
amount of modes. We describe the dynamics of this sys¬ 
tem by solving the Schrodinger equation 

H\ip(t)) (38) 

which can be set in matrix form as 


are not coupled to the excited state in the case where 
we consider modes with a same energy to. We are thus 
free to replace the whole set of modes by an effective 
state | A) characterized by an effective coupling constant 
which obeys 

Geff = v /|Gi| 2 + |G 2 | 2 + --- + |G jv | 2 , (34) 

that can be considered equivalent to the summation rule 
of the T factors over various decay channels which is tra¬ 
ditionally used in the context of the Fermi golden rule in 
QED and particle physics as well. However, in order to 


c e (t) 

d c gd (/) 
d t ■ 

_ C gWcomb W. 


-iw eg -G* (wi) ... -G* (wAT comb )' 


c e (t) 

G (u i) — iwi 

0 


Cg,l (t) 

-G(u ncomb ) 0 -\uj v comb . 


_ C g,iV comb (*). 


( 39 ) 
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l c g,A(0| 2 |c e (f)| 2 




Figure 1. Plot of |c g ,A (t)| 2 for various modes for a system 
with A^ C o m b = 100 resonant modes with frequencies within 
A to = 16 r of w eg . Note that at short times the behaviour 
is quadratic for all modes: this reflects the initial, coherent 
response of the electromagnetic vacuum. 


The eigenvalues (= k\) and eigenvectors (= |ka)) of H 
are estimated numerically and thereafter are used to solve 
for c e ( t ) and c gj A (t): 


c e ( t ) 

c g,i 0 ) 
Cg.Afcomb (0 


JV+1 

= X W Kj)e 

A—1 


(40) 


To describe a situation in which the atom is initially ex¬ 
cited, in a zero temperature environment, we consider 
the initial conditions 


/ c e (t = 0) = 1, 

\ Cg,A (t = 0) = 0 


(41) 


B. Numerical results 

We now consider -/V com b electromagnetic modes form¬ 
ing a frequency comb around w eg . The frequency spacing 
between two successive modes is taken to be much smal¬ 
ler than w eg . If one wants to generate the dynamics of 
the Fermi regime, one must make sure that the total fre¬ 
quency range Aw of these electromagnetic field modes 
satisfies 


Aw > r 


(42) 


with T given by (161, which can be seen as a manifesta¬ 


tion of the time-energy uncertainty AEAt ^ H/ 2. 

If Aw is taken to be much smaller than advised by 
then some modes which are significantly excited by their 
interaction with the atom are wrongly ignored. If Aw 
is taken to be much larger, then off-resonant modes are 
needlessly included in the treatment [TO]. One can indeed 
distinguish the “active” from the “inactive” modes by 
looking at the |c gj A (t) | 2 curves of Fig. [lj 


Figure 2. Transition towards exponential decay by fulfilling 
better and better the condition Aw T. 


In practice, we varied the comb width Aw in order to 
look for the change in behaviour when Aw approached 
a value close to F. As to the probability amplitudes 
|c gj A (t)| 2 of the field modes, one would expect from the 
expression for |c gi A (t )| (see (25)) to see a quadratic curve 
at short times. Notice that this initial behaviour should 
be entirely independent of frequency. These expectations 
are nicely met in Fig. [I] Fig. [2]shows how one gets closer 
to the Fermi regime -which constitutes the beginning of 
the exponential decay- by widening the electromagnetic 
frequency comb. 


We used, for our numerical calculations, the exper¬ 
imental values for the 2p — Is transition in atomic 
hydrogen, namely, w eg = 1.55 x 10 16 s -1 , |d ge | = 
6.31 x 10 -30 C m and l/ r = 1-60 x 10" 9 s [□]. The 
mode discretization described above and in the App. is 
quite robust in terms of its dependence on N c om b, as illus¬ 
trated by Fig. [3| where we compared the survival prob¬ 
ability |c e (t) | 2 found from numerical calculations to the 
standard result of Wigner-Weisskopf perturbation theory. 

The survival probability displays a “revival” due to 
the discretization of the outgoing modes, which results 
in a finite recurrence time t reC ur- However, a number 
of iV com b = 100 quasiresonant modes is enough to send 
this recurrence time to around 10/r (see Fig. [3]), which 
is more than enough for what we want to study here. 
Even with IV com b = 20, the discrepancy with the W-W 
theory in the region of interest (say 0<t<2.5/r« 
4 x ICE 8 s) is not noticeable. Note that the condition 
for the “weak coupling limit”, T <C w eg , is de facto met 
for the system we considered. We can also detect on 
Fig. i (see curve with Aw = 16 T) the initial Zeno regime 
after which the curve indeed approximates the Wigner- 
Weisskopf exponential. We shall consider later the role 
played by high-frequency, off-resonant modes in the Zeno 
zone. 

Our results agree with the Fermi golden rule in the 
interval w” 1 < t < F but also with the Wigner-Weisskopf 
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MO | 2 \CgAt )\ 2 




Figure 3. Decay of the survival probability |c e (f) 2 for calcu¬ 
lations with jV comb = 20 (dotted blue), N colnb = 40 (dot- 
dashed red) and iV com b = 100 (dashed green) against the 
theoretical expected behaviour exp (—Ft) (solid black). The 
curves for IV = 20, IV = 40 and IV = 100 are virtually indis¬ 
tinguishable for t < 2.5/r (remember 1 /V = 1.60 x 10~ 9 s). 


perturbation theory for longer times. For N comb = 100 
quasiresonant modes, we plotted in Figs.[4]and[5]the time 
evolution of |c gi A it) | 2 for two modes, one rather close 
to the resonant transition frequency (Fig. HI), and the 
other one extremely close to resonance (Fig. [5]). Wigner- 
Weisskopf theory predicts jT2] 


c g ,A (t) = e IWA G a (r 0 ) 


3 -i(w eg -w A )t e -rt/2 _ I 
Wa — W e g + ir/2 


(43) 


which tends to 


Figure 4. Evolution, for a system with iV comb = 100 
quasiresonant modes of the field, of the probability that the 
100 th mode of the held is occupied (solid red), plotted against 
the expected Wigner-Weisskopf behaviour (431 (dotted blue) 
and the asymptotic limit ( |44[ ) of the latter expression (dashed 
black). 


ICg.A(0 | 2 



]Cg,A (t) 


12 rt—»-+c 


|Ga (ro) 


(w A - w eg ) 2 + (r/ 2 ) 2 


(44) 


as t —> +oo. These expectations are nicely met too in 
Figs. [4] and [5] where the asymptotic limits (44) are expli¬ 
citly drawn. 


Figure 5. Evolution, for a system with iV comb = 100 
quasiresonant modes of the held, of the probability that the 
(very resonant) 50 th mode of the held is occupied (solid 
red), plot ted against the expected Wigner-Weiss kopf beha¬ 
viour ( |43| > (dotted blue) and the asymptotic limit (441 of the 
latter expression (dashed black). 


We investigate in several steps the influence of off- 
resonant electromagnetic modes on the decay of the sur¬ 
vival probability. 

First (Sect. IVI, we use the E • x interaction Hamilto¬ 
nian in the dipole approximation of Sect. m and dis¬ 
card all frequencies above the cutoff frequency cKq of 
Sect. m while in Sect. [V] we make use of the exact 
matrix elements for the E • x interaction Hamiltonian. 

Then (Sect. VI), we repeat the same proce ss with the 
A • p interaction Hamiltonian. In Sect. |VI A we use the 
dipole-approximated coupling derived similarly to that 
of Sect. IIA and discard all frequencies above the cutoff 
frequency cA'c of Sect. |II C[ while in Sect. |VIBl we make 
use of the exact matrix elements for the A ■ p interaction 
Hamiltonian found in 1131. 


IV. E x COUPLING IN THE DIPOLE 
APPROXIMATION 


We focus here on the dipole approximation. Remem¬ 
ber that the integral (25) diverges at all times, which 
requires us to introduce a cutoff frequency, as was done 
in Sect. IIC If we follow this prescription, the integral 
(25) becomes 


= -Wsmc 2 ((««-“) \ 

(45) 

This expression of the decay probability is compared with 
the result of the usual Wigner-Weisskopf approximation 
on Fig. [6} One notices, at very short times -the so-called 
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Figure 6. Linear-log plot of the survival probability 
|c e (t)| 2 from Wigner-Weisskopf theory (dashed black) and 
the truncated cardinal sine integral ( |45[ ) (dotted orange) with 
|Ce(t)| 2 = l — Pde cay (t) . The two main time regimes (strong 
emission at very short times and exponential decay later) 
can be visualized on this single graphic. Remember l/wc = 
2.09 x 10“ 2O s, 1/weg = 6.45 x 10“ 17 s, l/V = 1.60 x 10“ 9 s. 


Zeno regime- a strong drop during the quadratic decay of 
the truncated cardinal sine integral (451. At longer times, 
the integral (45) follows the linear decay predicted by 
Fermi’s golden rule -with an offset- and fails to reproduce 
the predictions of the Wigner-Weisskopf approximation 
beyond t > 1/F, which is expected from this first-order 
approach. As we shall show this discrepancy can be sup¬ 
pressed by going beyond first order perturbation theory 
to describe quasiresonant electromagnetic modes. 


Our goal here is to connect the two regimes (Zeno for 
short times and W-W for long times) observed on Fig. [6] 
For intermediate times, where the Fermi golden rule is 
valid, we know that first order perturbation prevails. Be¬ 
sides, it is notoriously known that perturbation theory 
is not accurate when we consider non-perturbed states 
with very close energies. In order to develop a perturb¬ 
ative approach, we treated the quasiresonant modes first 


by diagonalising the Hamiltonian as done in Sect. Ill 


before incorporating the off-resonant modes of the elec¬ 
tromagnetic field through time-independent perturbation 
theory. Let us denote by Ho the Hamiltonian which de¬ 
scribes the atom, the field, and the interaction between 
the atom and the quasiresonant modes of the electromag¬ 
netic field, and Hhigh the Hamiltonian which describes 
the interaction between the atom and the other (off- 


resonant) modes. Namely 


ATcomb “I - Afoffshell 

H 0 = fiw e g|e,0)(e,0| + 55 huj x \g, 1a)( g, 1*1 

A—1 

Ncomb 

+ i^ 55 [Ga ( r o) |g, !A>(e, 0| — GX (r 0 ) |e, 0)<g, 1 A |] 

A—1 

(46a) 

and 


ATcomb “1“ AToffshell 

H h igh = ifi 55 [ G A(r 0 )|g, lA>(e,0| 

A —iVcomb “1“ 1 

( r o) |e,0)(g,l A |] 


(46b) 


where lV C omb is, as previously, the number of quasireson¬ 
ant modes and -/V of j she n is the number of off-resonant 
modes. The effective frequencies uj\ and couplings 
G A (ro) are defined in App. 2 a We wrote |g, 1a) = 
|g) <S> | 1a) where |g) is the atomic ground state while the 
single photon-single mode states |l A ) ar e obt ained by the 

action of the creation operators (havg)^ (A.9) on the elec¬ 
tromagnetic vacuum |0). Now, assume that we diagon¬ 
alized Hq- Since the off-resonant modes are uncoupled, 
the spectral decomposition of H 0 reads 


A^comb AT CO mb “h-A^offshell 

h 0 = 55 iV’iXV'ii + Hujx | g, 1a) <g, 1a | 

2=0 A —iVcomb+ 1 

(47a) 

with 


JV C , 


JV C , 


IV 7 *) = e o Ie, 0) + 55 e A I S> x a) = 55 c aI a )> ( 47b ) 


A=1 


A=0 


where we introduced the new notatio ns |e , 0) = 10), and 
|g, 1a) = |A). It is useful to notice that ( 47b| ) can be inver¬ 
ted by 


A^comb 


l A ) = 55 ( e i) l^>- 


A—0 


(48) 


The e\ coefficients in (47b) were obtained numerically in 
the previous Sect. Ill by diagonalising the matrix (28). 
Expanding the relations (A | £) = S\^ and (ipk \ ipj) = Sjk 
yields, respectively, the useful relations 


ATcomb 

f: ( c a) 6 c = 

i =0 

A^comb 

E W)‘4 = v 

A—0 


(49a) 

(49b) 


Treating Hhigh to first order in perturbation we get the 
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normalized perturbed eigenstates 

ATcomb “1“AToffshell 


l^} (1) = 


1 


VNl 


Ik i\) (1) = 1 

IS ’ A ' 

with 


i / \ . • i (r 0 ) 

1^) + 1£ ° E ( n i_ n \ 

A=JV comb + l K W a) 


— l + k? 


Is,1a > + i (Z7^kfr |e, °) 

V^eg j 


A^comb AT 0 ffshell I / \ 1 2 

|2 Y" |Ga (rp)| 

1 ^ / a , \2’ 

A=JV comb ( w 0 ~ w Aj 


Is, !a) 

(50a) 

(50b) 


JV A = 1 + 


|Ga (ro)|l 

(w eg -W A )' 


(51a) 

(51b) 


Note that to first order in G\ (r 0 ) / [(w eg — wa)], these 
perturbed states are orthogonal. This means that, still 
to first order in these parameters, the perturbed spectral 
decomposition of H = H 0 + Hhigh and the corresponding 
evolution operator are valid. 

The initially excited state | tp [t = 0)) = | e, 0) can be 
expanded over the perturbed eigenstates, making use of 
d49b and (pi): 


w ‘ = °»-7? 


Afcomb 

E Ev*(4)*l^> (1) 

L A=0 

ATcomb “1“ AToffshell 

E 

A — -ATcomb - ! - 1 


Ga (tq) 

(Weg - W A ) 


|g, 1 a> (1) 


1 

Vn 


A^comb“1“A^offshell I / \|2 

1+ |Ga ( r o)l 


A—Af com b +1 


(cJ e g Cd\) 


and thus 


A^comb^AToffshell \ ( \ 12 

cjv = i+ y |G>(r ° )l 

(w eg -u>x) 


|e,0) 

(52) 

(53) 


A—AT comb + l 


Then we can compute the survival probability at time 
t. For that we let the evolution operator U ( t ) act on 
| <p (t = 0)), which yields 


\<P(t)) = 


Vn 


ATcomb 


2 — 0 


A^comb “1“ A^offshell 

-i E \/E 

A—A/comb +1 


S ,1a) (1) 


_Ga(to) iWAt 


(w e g UJ\) 


(54) 


The survival probability is then deduced with the help of 
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Figure 7. Decay of the survival probability |c e (f)| 2 with the 
dipole-approximated E ■ x coupling from numerical calcula¬ 
tions (solid blue), Wigner-Weisskopf theory (dashed black) 
and the truncated cardinal sine integral ( |45| ) (dotted or¬ 
ange) with |c e (f)| 2 = 1 — Pdecay (f). The logarithmic scale 
for time is used so that the two main time regimes (strong 
emission at very short times and exponential decay later), 
as well as the transition between them, can be visualized 
on this single graphic. Quasiresonant modes are discret¬ 
ized with N com h = 100 and off-resonant modes are discret¬ 
ized with IVoffsheii = 100. Remember l/u>c = 2.09 x 10“ 20 s, 
l/o; eg = 6.45 x 10 _17 s and 1/F = 1.60 x 10~ 9 s. 


(50): 


-Psurv (t) = |(e, 0 | </? (t))\ 


1 

N 


N Ci 


E 

2 — 0 




A^comb “1“ A^ 0 ff s hell I / \ I 2 

A = JVcom b + l ^ e S “ 

(55) 

To compute this quantity we resort to numerical calcula¬ 
tions. The main goal here is to link the initial behaviour 
of the survival probability as given by the cardinal sine 
integral (25) to the experimentally observed exponential 


decay at large times. The expression (55) of the survival 
probability is compared in Fig. [7] with the usual Wigner- 
Weisskopf exponential exp (— Ft) and the truncated car¬ 
dinal sine integral (|45|). 


We see on Fig. [7] that our method allows the numerical 
results to espouse, at short times, the predictions of first- 
order perturbation theory, and, at large times, the usual 
Wigner-Weisskopf exponential. As we shall show in the 
next Sect. [Vj however, the fast drop at short times seen 
in Figs. [6] and [7] is not physical. 
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V. EXACT E X COUPLING 
A. Results 

When going beyond the dipole approximation, one can 
derive a more exact expression for the matrix element: 


GA(ro) = iV^ dse ' vr(ro) 


arctan 


“o + 


(a)' 


1+ ^ 


(a)' 


1+ ^ 


(a)' 


WA 

UJx 


(56) 


with natural cutoff frequency wx = (3c) / (2ao), a result 
which, as far as we know, is original. This is done in the 
next Sect. | V B| This time, first-order perturbation theory 
yields the following cardinal sine integral 


(f) = l d “ " s ”“ 2 ( ( "«- 0 

-i 2 


(a)' 


1+ ;^ 


3 

2 + 2 


arctan ( 


(a) 


1 + 7 ^- 


(a)' 


(57) 


i/ojx Loj e g i/r 



Figure 8. Decay of the survival probability |c e (t)| with the 
exact E ■ x coupling from numerical calculations (solid blue), 
Wiener- Wei sskopf theory (dashed black) and the cardinal sine 
integral (571 (dotted orange) with |c e (f)| 2 = 1 — Pdecay (t). At 
long times numerical calculations are indistinguishable from 
Wigner-Weisskopf decay. A logarithmic scale for time is used. 
Quasiresonant modes are discretized with N com b = 100 and 
off-resonant modes are discretized with A 0 g s h e ii = 100. Re¬ 
member l/wx = 1.18 x 10 _19 s, 1/cUeg = 6.45 x 10~ 17 s and 
1/r = 1.60 x 10 -9 s. 


canonically conjugate momentum to the vector potential 
is given by 


The same perturbation method as that exposed in the 
previous Sect. |IV| is used. The frequencies oj\ and coup¬ 


lings G\ (r 0 ) are defined in App. 2 b Results with this 
coupling are shown in Fig. [8j 

As noted above, the fast drop in the survival probabil¬ 
ity is an artifact of the truncated dipole approximation. 
Note that our perturbation method again gives precise 
results at both short, intermediate, and long times (see 
Fig. [8}. 


2 

n(x) = iE 

x=l 


d 3 k I he |k| 

(27t) 3 y 2eo 


(&(*) (k) e lk ' r e (>c) (k) 



Finally the interaction Hamiltonian can be cast in the 
following form: 


B. Derivation of the interaction matrix element 


We now prove (56). In the Power-Zienau-Woolley 


picture of quantum electrodynamics, the interaction 
Hamiltonian reads m 


Hi - 


1 

eo 


dxP(x) -n(x). 


(58) 


The electric polarization field reads, in the case of atomic 
hydrogen 


P (x) = —er / du<5(x — ur) (59) 

Jo 


where r is the electron position operator and e the ele¬ 
mentary electric charge (charge of a proton), and the 


Hi = i e 



+oo 

«(*) ( k ) e (x) ( k ) 

n—0 


(ik-r)" 

(n + 1)! 



(61) 


To compute the matrix element of this interaction 
Hamiltonian between the Is and 2p levels, the hard 
part is to compute (Is | (ik-r)”r | 2p77i 2 ) where m 2 
is the magnetic quantum number for the 2p sublevel 
considered. We follow the first steps given in [13] for 
the derivation of the corresponding interaction matrix 
element in the A ■ p coupling. 

Remember the expressions for the wave functions of 












































11 


the Is and 2 pm 2 sublevels: 


Ipls ( x ) = 


exp 


(- 2 ) 


V’2p m 2 ( x ) — 


ex p(-S) 72 


' ttOq «o 


x ' £m 2 • 


The t m2 are given by 


£o = e z, 

o, rlz ie„ 
£±i = =F- 


72 


The three states (62b) span the 2p subspace for the elec¬ 


tron, such that any state on the 2 p sublevel is a linear 
superposition of these three states. Choosing a coordin- 
(62a) ate system for which k/ |k| is the third basis vector, we 
compute 


(62b) 


(63a) 

(63b) 


(ls| (ik- r)”r| 2 pm 2 ) = j dx (ik • x) n x (x) (x) 

- / d<p f AO sin<9 cos n 9 [ + dx exp —) x n (x-£ TO2 ) x. 

Jo Jo Jo \ 2 oq / 


(64) 


87 t a o 


In our basis where k/ |k| is the third basis vector, we have 



sin 6 cos ip 

(^m 2 sin 9 cos ip + ^ sin 9 sin p + cos | x | 2 

sin 9 sin p 
cos 9 


(65) 


The stand for the components of £ TO2 in that basis. After carrying out the easy integrations over ip and changing 
variables (rj = cos 0 ), we get 


(Is | (ik • f)" r 12 p m 2 ) = 77 =^ [ <7 
4y2aQ J -1 


S™-2 V I J 

(ml ( 1 - 7 ) 

2^71 7 


da; a 4+ " exp - 




2 ao 


( 66 ) 


Since, according to (61), the dot-product of ( 66 ) with the polarization vectors e.t x \ (k) is to be taken, its component 
along k/ |k| can be discarded !}I3]. One can thus simply compute the scalar quantity 


(ls| (ik- r) n r|2pm 2 )_L = (k) [ dr] (l - 7) 17 ” [ dxa 4+n exp ) 

4v2a 0 J —1 J 0 \ " ^0 / 


a o 
1 1 
272 a 0 


l + (-l) r 


■ 4 ^ ^ (1 + n) (3 + n) ^2a 0 


g \ — (n+5) 

(n + 4)! £ m2 • e (jf) (k). 


Then we compute the sum over all values of n, according to the expression found in (61): 


Is 


00 ... n 

y ( lk ■ r ) j. 

n 7+1)! 

n—0 v 7 


9 \ _ $m 2 • e (*r) (k) / 3 \ y— ,. ,, isn 1 + (~ 1 ) ( 3 ^ {n + 4 )! 

/ ± 272<7 \2aoJ ^ (l + n)(3 + n) \2a 0 ) 7 + 1)! 

=(3 ) 5 - •*■> « t (3 ' k ' *)* (- 1 )" ■ 


(67) 


( 68 ) 


It then comes in handy to notice that 

7 4 - 1 ) (71 + 2 ) 
( 2 n+l) 


(n + 1 ) + - ( 1 + 


1 


2 n + 1 


(69) 
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to rewrite (681 as 


Is 


00 /‘I ''\ n 

v (lk • r ) r 

n (n+ 1)! 
n=0 x 7 


2p m 2 ) =\/2(-j a 0 £ m2 ■ e M (k) 


!+ (|a 0 |k|) 



■ (70) 


Using (611 and the following relations m 


1 110c 


2^ 

dge = 6 &0 gg £m 2 ; 

_ 3 h 1 

Weg ~ » 7W ^2 

o 777/g a 0 


(71a) 

(71b) 


we then retrieve (561. 


VI. A p COUPLING 

A. Dipole approximation 

In the dipole approximation, the matrix element of the 
interaction Hamiltonian A ■ p is m 


G\ (r 0 ) — 


to, 


eg 




dgo • vf (r 0 ). 


(72) 


This time, first-order perturbation theory yields the fol¬ 
lowing cardinal sine integral 


P, 


decay 


Id 

(t)= 1 s ° 


(O a 


6ir 2 eohc 3 eg 
cK c 


r cK c / 

J d oj uj sine 2 ^(w eg — to) -J . (73) 

The same perturbation method as that exposed in 
Sect. |IV| is used. The fr eque ncies u\ and couplings 
G\ (r 0 ) are defined in App. |2c| Results with this coup¬ 
ling are shown in Figs. [9] (short times) and [To] (long 
times). A comparison of Figs. [9] and [7] shows that, even 
in the dipole approximation, the Zeno effect obtained is 
very much reduced in the A ■ p coupling compared to 
what it is in the E • x coupling. 


B. Exact A ■ p coupling 

When going beyond the dipole approximation, one can 
derive a more exact expression for the matrix element, 
namely jl3| 


G x (r 0 ) = 


u eg 


x/2he 0 u\ 


dge • Vp 


( r o) 


1 


(s)' 


1 + ^ 


(74) 



Figure 9. Decay of the survival probability |c e (t)| 2 with the 
dipole-approximated A ■ p coupling from numerical calcula¬ 
tions (solid blue), Wigner-Weisskopf theory (dashed black) 
and the truncated cardinal sine integral ( |73| ) (dotted orange) 
with |c e (t)| 2 = 1 — Pdecay (t). A logarithmic scale for time is 
used. Quasiresonant modes are discretized with N coln b = 100 
and off-resonant modes are discretized with A 0 fj s h e ii = 100. 
Remember l/wc = 2.09 x 10~ 20 s, l/io eg = 6.45 x 10~ 17 s 
and 1/r = 1.60 x 10 _9 s. 

with natural cutoff frequency wx = (3c) / (2ao). A mode 
discretization similar to that used in Sect. [iVj and the 
important results given in App. 1 2 d[ This time, first- 
order perturbation theory yields the following cardinal 
sine integral 


t 2 Id. 


gel 


T^decay 

r°° 

/ dw 


1 + ^ 


(*)' 


1 4 


sine 


(w eg - to) - . (75) 


The same perturbation method as that exposed in 
Sect. IV is used, with the following modifications. The 
frequencies oj\ and couplings G\ (ro) are defined in 
App. 1 2 d| Results with this coupling are shown in 
Figs. [IT (short times) and [l2] (long times). Figs. [8] 
through in show that the strong decay of the survival 
probability observed on Fig. [7] for very short times is an 
incorrect prediction. This highlights the fact that, if the 
dipole approximation is to be used at all, this should be 
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1 lUJQ 1/w e g 1/f 



Figure 10. Decay of the survival probability |c e (t )| 2 with the 
dipole-approximated A ■ p coupling from numerical calcula¬ 
tions (solid blue), Wigner-Weisskopf theory (dashed black) 
and the truncated cardinal sine integral ( |73| ) (dotted orange) 
with | c e (£) | 2 = 1 — -Pdecay (t) . At long times numerical cal¬ 
culations are indistinguishable from Wigner-Weisskopf de¬ 
cay. A logarithmic scale for time is used. Quasiresonant 
modes are discretized with lV com b = 100 and off-resonant 
modes are discretized with N 0 f[ s h e ii = 100. Remember 
1/wc = 2.09 x 10 _2 ° s, 1/weg = 6.45 x 10 _17 s and 1/F = 
1.60 x 10“ 9 s. 


done with much care. Note that for short times our per¬ 
turbative approach fits the analytical results derived in 
m, by a “tour de force” which allowed the authors of 
m to estimate departures from exponential decay for 
short but also very long times in an exact fashion. 

It is worth noting that while Facchi and Pascazio con¬ 
sidered m the ratio between the Zeno time rz and 
the lifetime te of the excited level, where P surv (t) ~ 

t-t o 

1 — (t/r z) 2 and = 1/r, as the relevant parameter con¬ 
cerning the experimental observability of the Zeno re¬ 
gion, we argue that the relevant ratio is that between 
the “cutoff time” ix and the Zeno time rz. The cutoff 
time is simply defined as rx = 1/wx with wx the cutoff 
frequency of the atom-field interaction. Indeed, the re¬ 
sponse of the electromagnetic modes is no longer coherent 
after a time of order rx (as can be seen for instance on 
Fig. § where the cutoff frequency is wg). This means 
that after rx, one exits the Zeno regime in which the 
survival probability decays quadratically. With this in 
mind, we take an interest in the scaling properties of the 
ratio rx/rz- For hydrogen-like atoms with Z protons, 
the Zeno time scales, as noted in [15], like Z~ 2 . For such 
systems, since the Bohr radius scales like Z _1 , the cutoff 
frequency scales like Z , and the cutoff time like Z _1 , 
so that the ratio rx/rz scales like Z, while Facchi and 
Pascasio’s parameter (■tz/te) scales like Z 2 . 


Muj x 


1/OJgg 



Figure 11. Decay of the survival probability |c e (f)| 2 with the 
exact A ■ p coupling from numerical calculations (solid blue), 
Wigner- Weis skopf theory (dashed black) and the cardinal sine 
integral (751 (dotted orange) with |c e (t)| 2 = 1 — Pdecay (t). 
Numerical calculations are indistinguishable from (751. A log¬ 
arithmic scale for time is used. Quasiresonant modes are dis¬ 
cretized with Aomb = 100 and off-resonant modes are discret¬ 
ized with IVoffsheii = 100. Remember 1/wx = 1.18 x 10~ 19 s, 
1/weg = 6.45 x 10~ 17 s and 1/r = 1.60 x 10 _9 s. 


VII. CONCLUSION AND OPEN QUESTIONS 

We scrutinized the onset of the exponential decay of 
a quantum dipole in the Zeno regime. We focused on 
the role played by excitations of the off-shell electro¬ 
magnetic modes that are usually neglected. Our main 
result is that, before an exponential behaviour is mani¬ 
fest and even before the linear in time decrease of the sur¬ 
vival probability characteristic of the Fermi golden rule 
appears, a sudden decrease occurs for very short times, 
quadratic in time, which is directly related to the so- 
called quantum Zeno effect. 

As can be seen in Fig. [l3] (where the onshell and offshell 
contributions to the decay of the survival probability as 
expressed by the integral ( [57] ) are plotted separately), 
the offshell modes play a crucial role in the Zeno regime. 
They are activated very quickly but after a time of the 
order of the inverse of the cut-off frequency they do no 
longer respond coherently and their contribution to the 
decay vanishes. The linear behavior predicted by the 
Fermi golden rule occurs much later and is due to onshell 
contributions only. 

As a by-product of our analysis, we developed a per¬ 
turbative scheme that interpolates between the Zeno 
(parabolic) regime and the W-W (exponential) regime, 
passing through the Fermi (linear) regime. In the case 
of the 2p — Is transition we cover a very large interval 
of time (from 10“ 19 to 10 -9 second). Our results fit 
with the exact solution of Facchi and Pascasio m if we 
treat the A • p coupling without approximation, but our 
perturbative approach also makes it possible to treat the 
exact E ■ x coupling as well as the A ■ p and E ■ x coupling 
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Figure 12. Decay of the survival probability |c e (t)| 2 with the 
exact A ■ p coupling from numerical calculations (solid blue), 
Wiener- Wei sskopf theory (dashed black) and the cardinal sine 
integral (751 (dotted orange) with |c e (t)| 2 = 1 — Fdecay (t). At 
long times numerical calculations are indistinguishable from 
Wigner-Weisskopf decay. A logarithmic scale for time is used. 
Quasiresonant modes are discretized with A com b = 100 and 
off-resonant modes are discretized with A 0 ft s h e ii = 100. Re¬ 
member 1 /wx = 1.18 x 10~ 19 s, l/weg = 6.45 x 10 _17 s and 
1/F = 1.60 x 10 -9 s. 


1 lUJX tWgg i/r 



lime(s) 


Figure 13. Decay of the survival probability |c e (t)| 2 for the 
exact E-x coupling if one takes into account either only reson¬ 
ant electromagnetic modes (dashed green), only off-resonant 
modes (solid violet) and all modes (dotted orange). Re¬ 
member 1 /wx = 1.18 x 10 _19 s, l/weg = 6.45 x 10 _17 s and 
1/F = 1.60 x 10 -9 s. 


in the dipolar approximation. 

It is natural to ask whether the Zeno effect predicted in 
our paper is likely to be observed experimentally. The ef¬ 
fect that we predict, in the case of the ls-2p transition is, 
indeed, very small and does not occur during a measur¬ 
able time period (typically, of the order of an attosecond). 
This explains why our counter-intuitive predictions have 


not yet been verified experimentally. This verification, 
however, is not impossible in principle, but requires 

1 . to measure the change of population of the emitters 
carefully. 


2 . to investigate other unstable quantum systems 
(magnetic dipoles [9, IT?] , solid state or NMR 
qubits, Rydberg states [18] and so on). The reason 
is that the duration and intensity of the Zeno ef¬ 
fect is strongly dependent on the off-shell coupling 
parameters while the W-W and Fermi golden rule 
regimes are dominated by the on-shell coupling. Up 
to now, we limited ourselves to the study of a single 
transition (the Is — 2p transition for atomic hy¬ 
drogen) but we are convinced that there exist in 
Nature other transitions (quantum jumps) where 
the particular type of Zeno effect predicted by us 
(collective response of the off-shell modes at zero 
temperature) is enhanced. Magnetic dipolar trans¬ 
itions 0111] for instance are characterized by very 
long life times (in some cases a million times longer 
than the Is—2p transition), which suggests that the 
Zeno time could approach the picosecond, a tem¬ 
poral regime which can be probed by sophisticated 
pump-probe techniques such as the so-called ul¬ 
trafast transient absorption spectroscopy in which 
very quick responses (of the order of 200 fs) are 
reachable experimentally. 


Another interesting question concerns causality. It is 
known that, in accordance with Hegerfeldt’s theorem, 
even if we stick to the W-W approximation, a non-zero 
component of the single photon wave function will be 
present outside from the light cone expanding centered 
in spacetime on the atom at time t = 0 jT9ff2Tj . The 
reason therefor is that only positive spatial frequencies 
are present in the spectrum of this wave function, so that, 
in virtue of Paley-Wiener theorem, the domain of the 
wave function cannot be bounded. We investigated the 
departures from causality in the W-W regime in another 
paper and found that they are extremely small |22[. 

One could wonder whether more significant violations 
of Einsteinian causality might be present in the Zeno re¬ 
gime. In order to tackle this question, it is necessary 
to have at one’s disposal extremely accurate estimates of 
the wave function of the emitted photon. This is possible 
in our approach which provides an accurate estimate of 
the solution c e (f) of the equations of motion (141 that in¬ 
terpolates between the Zeno time and the lifetime of the 
dipole. This suffices in order to compute the f ield e mitted 
by the atom. Indeed, the formal solution of (14b) reads 


Cg,A (t) = Gx (r 0 ) f Ce(t') e^-^' At'. 
Jo 


(76) 


The information about the single-photon field (and in 
particular of its Glauber wave function) emitted by the 
atom is entirely encoded in the coefficients c gj >, (t). We 
expect therefore to be able to study in depth possible 
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violations of Einsteinian causality for very short times. 
Of course, we performed several approximations in our 
approach (RWA for instance) and it remains an open 
question to know whether violations of causality are due 
to the approximations - as suggested in 0 - or whether 
they are unavoidable features of QED. These questions 
are complex but at least our approach makes it possible 
to tackle the problem in an original manner. 


Appendix: Mode discretization: energy slice ladder 
operators 

1. Position of the problem 


In the present appendix we show how to discretize the 
continuous modes of free space through the introduction 
of creation and annihilation photon operators for energy 
slices. We are interested in the spontaneous emission 
of light in a vacuum and hence our starting point will 
feature a continuous infinity of modes. Hence we first 
want to adapt the equations of Sect. |II A| to this situation. 
In Sect. |II C] we gave part of the prescription for switching 
from a discrete sum over arbitrary modes to a continuous 
sum over the free space modes. This was done by writing 


1 

V 


E-E 


d 3 k 


(A.l) 


where >r are the two possible photon polarizations, V is 
the initial quantization volume and the modal functions 
are given, in the Coulomb gauge, by 

v a ( r ) (k, r) = e M (k) . (A. 2) 


It will come in handy here to choose linear polarization 
vectors as a basis. These are given by 


6 (q) (k) = cos^cos p e x + cos(9sin<p e y — sin0 e z , 
6 ( 2 ) (k) = — sin p e x + cos p e y 

(A.3) 

where 9 and </> are the spherical coordinates of the wave 
vector k with respect to an arbitrary cartesian frame 
(here, we will take e~ in the direction of d ge . A pres- 
ciption for going over to the continuous case for ladder 
operators needs to be added. It consists in the substitu¬ 
tion 


a A) a^, =6 AA /-s> a* (k), (q) = (2 tt) 3 5 (k - q) 5^. 

(A-4) 

Now we want to be able to perform numerical compu¬ 
tations, and hence have to discretize this continuum in 
some way. 

2. Results 


a. Detailed treatment in the case of the dipolar E ■ x 
coupling 


We then rewrite the rotating interaction Hamiltonian 
( 9| w ith the use of (A.l I, (A.2) and, in the second step, 
( A.3| as 


H? = - —E E / d W d $ sind J — / d k hi [a x (k, 9, p) e (x) (9, <p) ■ d*Je)(g| - H.c.] 

(27T) x=1 Jo "'° » ze ° Jo 


(2ir Y 


^ |d -' 


dp / d0 sin0 / dk k 2 (— sin0) ai(k,6,(p) \e)(g\ —a\(k : 6,ip) \g)(e\ (A.5) 


where we took d ge as the reference axis for the polar coordinate 9, and used the fact that, for the 2 p — Is transition 
in atomic hydrogen, the dipole matrix element is real m- With the prescription to obtain an effective coupling given 
below (321 in mind, we introduce the spherically averaged creation operators 


p2n pn / p2n p tt 

/ dp d 9 sin 2 9 d\ (k, 9, p) = a sp h (k) ( dp d 9 sin 3 9 

Jo Jo \J o Jo 


2n 

= 2 y “^“dsph (k). 


(A.6) 


r 


In the leftmost member of (A.6), sin“ 9 is the product of 
a first sine coming from the differential solid angle ele¬ 
ment, with another one coming from the coupling (see 
(A.5). According to our prescription, the coupling of one 


should be squared to obtain, after summing over a given 
ensemble of modes and taking the square root, the effect¬ 
ive coupling constant of that ensemble, hence the third 
power of sin0 in the central term of (A.6). With these 


particular mode (here, with angular orientation [9,p)) new operators, the interaction Hamiltonian (A.51 can be 
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rewritten 


H? = 


2i 


(2^r 


hen 


3en 


5 

dk k 2 


«s P h (k) |e)(g | —«g ph (k) | g)<e | . (A.7) 


following commutation relation: 

[aavg>CU = <5ac- 


(A.11) 


From this we can compute the matrix elements of the 
interaction Hamiltonian between the excited state | e, 0) 

and the energy slice field state (davg)^ I g, 0), which reads 


It can be checked from (A.4) and (A .61 that the spheric¬ 


ally integrated ladder operators obey the following com¬ 
mutation relation: 


a sph (k), a f sph ( q ) 


(2nf 

kq 


6{k-q). 


(A.8) 


(e,0| H? (dav g ) f |g, °) = l d s e l { k x+i ^k\) 2 ■ 

(A.12) 

For a complete description we also need to compute the 
matrix elements of the free field Hamiltonian for these en¬ 
ergy slice field states. The electromagnetic field Hamilto¬ 
nian ([3| is rewritten, in free space, 


We further define 


i (27 t) 2 / dk k 2 a sp h (k) = a 

Jk x 


( r k \+i 

- y g ( J k dkk 3 


£.4 _ i- 4 ' 

K \+l K \ \ »A 
4 I “avg 


to rewrite (A.7) as 


(A.9) 


Hp = 


— IdgelV 

(2 7t) 2 V 3e 0 ' ge| y 


k 4 _ 1.4 \ 2 

E\+l ^A 


«avg I e ) (g I + (^avg |g)(e| (A. 10) 


where a* vg = Savg A+1 defined through (A.9) and the in¬ 
tervals [fcAj^A+i] are taken to have empty intersections 


and [0, +oo[ as their union. It can be checked from (A.8) 
and (A.9) that the energy slice ladder operators obey the 


H r = 


he 


E 


dep / dO sin 9 

Jo 


(2 kT^Jo 

pOO 

/ dfc A: 3 (k, 9, ip) (k, d, ip). (A.13) 

Jo 


As seen from (A.3), electromagnetic modes polarized 
orthogonally to the atomic dipole vector d ge are not 
coupled to the atom and thus evolve freely. It is thus 
only useful to consider 


(g>0| a avg & R aH |g, 0). 


(A. 14) 


*avg ~ 2 R a avg 

To that end it is best to first compute the commutator 


di (q, 0, ip) , (Savg) = - (2 tt) 


(K+i-KY 

sin 9 ^fq [0 (q - k\) - 0 (g - &a+i)] (A.15) 

where k\ and k\ + i are the bounds of the i th in terva l 
(ener gy sli ce). This can be seen to be true from (A.4), 
(A.6) and (|A.9). Noticing that since k\ + i > k \, we have 


[0 (q -k\)-@(q- fc A +i)] =Q(q- k\) + 0 (q - k\+i) - 20 (g - k\) 0 (g - k\+i) 


= © (g - k\) + 0 (g - fcA+i) - 20 (g - k\ +1 ) = 0 (g - k\) - 0 (g - k x +i ), (A.16) 


we compute (|A.14|) as follows: 

fg,0 


avg h r as. 


g ,°)=^«Ac 


3hc 

= ~2w ** 


= 4 he 5\£ 


f>2ir f-TT 


dip / dd sim 0 


d k -j-—- r (0 (g -k x )-Q(q- &a+i)) 

E\+l ^A 
rk a+i jA 


dk 


!k\ 


£.4 — k 4 

K A+1 K X 


dk 


!k\ 


k 4 — k 4 ^ 
E\+l ^A 0 


= 4 Sci 


k 4 — k 4 ’ 
^A+l ^A 


(A.17) 


r 


Notice that if one writes /ca+i — k x = Ak x , and lets The matrix elements (A.12) and (A.17) are the ones we 


Ak\ —► 0, then (A.17) tends towards hck\ ~ hck x+ \. 
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use for our numerical calculations. In the body of the 
article we denote them as (e, 0 | Hf- (dL g )^ | g, 0) = 
—i HG* X (r 0 ) for the coupling matrix (see (J9|) in the text) 


element and (g, 0 | a^ vg H R (a* vg ) f I S>°) = kwx for the 


avg“« Vavg 

diagonal field energy matrix elements. 


b. Results for the exact E ■ x coupling 


For the exact E • x coupling (74) we define 


F(k)= k ^x 


l+£ 


(*y 


arctan A- 


(A) 


1 + 


(A)' 


k 

kx. 


(A.18) 


and obtain 


<e,0|tff (a* vg ) f 


where 


M> = ^ yj I d ge I Weg 


(/W (fc A+1 ) -/ (4) (*a)) 5 (A.19) 


/ (4) (fc A+1 ) - j( 4 > (k x ) = / d k k 2 F 2 (k ) (A.20) 

Jk x 


with a^vg = a a v g A+1 defined through 


.(2tt)~ 


rk \+1 


d k k 2 F (fc) a sp h (fc) 


avg 




dfc k 2 F 2 (k) 


(A.21) 


c. Results for the dipolar A ■ p coupling 


For the A ■ p coupling in the dipole approximation (72) 
we obtain 


(e,0| H? (flavg) 1 |g,0) = ^y^|d ge | Weg 

h 2 — h 2 
^A+l ^A 


(A.24) 


with a* = a a v g A+1 defined through 


i (2 tt) 
and 


dk k 2 a sp h (k) = a: 


- a a 

avg 


dJfc kj 
(A.25) 


(§i 0|®a Vg Hr fiSglfcO) = 5 Af ^2 +1 _ ^2 • ( A - 26 ) 

o k a+i k a 


The matrix elements (A.24) and ( |A.26| ) are the ones 
we use for our numerical calculations. In the body of 
the article we denote them as (e,0\Hf (a* V g) |g, 0) = 
—iHG* x (ro) for the coupling matrix element and 

(g, 0|aav g Hr (da Vg ) f |g, 0) = hu}\ for the diagonal field 
energy matrix elements. 


d. Results for the exact A ■ p coupling 


Finally for the exact A • p coupling (741 we obtain 


gel “■'eg 


(e,0| H? (havg) 1 |g, 0) = ^ |d 

(/( 2 ) (k x+1 ) - / (2) (k x )) 


(A.27) 


We also find that 


(g, 01 a* vg Hr a^ g | g, 0) = He <5 Af 


where 


I (5) (fc A+ i) - J< 5 > (kx) 
'/W (fc A+1 )-/W (kx) 
(A.22) 


I (5) (fc A +i) - I (5) (kx) = f kX+1 dk k 3 F 2 (k) . (A.23) 

Jk x 


As far as we know the integrals (A.201 and (A.23) 


cannot be solved by hand and we compute them nu¬ 
merically. The matrix elements (A.19) and (A.22) 


are the ones we use for our numerical calculations. 
In the body of the article we denote them as 
(e,0|iJf (a^g) 1 |g,0) = —iKG* x (r 0 ) for the coupling 

matrix element and (g, 0|a* vg Hr (a* v g ) f |g, 0) = Sw A for 
the diagonal held energy matrix elements. 


where 


k 2 1 

lV)(k) = ^- 


6 


1+hSr 


(A)' 


n 3 


with fia Vg = aav g A+1 dehned through 

ki 


A avg — “avg 

P^A + l 


i (2 tt)' 


dk 


1+ A- 


(A)' 


^sph (&) 


^avg 


dk 




1+ hr: 


(A)' 


1 4 


(A.28) 


(A.29) 
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We also find 


where 


(g, 0|a* vg H r a^ g |g, 0) = he 5 XC 


/ 1 2 (3) 4 5 6 7 8 (k x+1 )-I^ (k x ) 
m (k x+1 ) - /(2) (k x ) 
(A.30) 


/ (3) (* a + i ) - I {3) (kx) = 


d k- 


k 2 


' k\ 


1 + £ 


(A) 


&x k X kx {kx ~ 3 k\) (k\ + 3 k x+1 kx {k\ - 3 k 2 x+1 ) (k' 
48 


Ui 


+ 3fcx) 


(A.31) 


—3arctan ( 

kx 


{kl + klY 


3 arctan 


C k 2 


A+l 


fc l) 


vA+1 

kx 


The matrix elements (A.27) and ( |A.30| ) are the ones 
we use for our numerical calculations. In the body of 
the article we denote them as (e, 0 |Hf (davg)^ |g, 0) = 
—i hG* x (ro) for the coupling martix element and 

(g, 0|da Vg Hr (a* vg ) f |g, 0) = hw a for the diagonal field 
energy matrix elements. 
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